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The cosmological backreaction from perturbations is clearly gauge-dependent, and obviously de- 
pends on the choice of averaged Hubble rate. We consider two common choices of Hubble rate and 
advocate the use of comoving volume-preserving gauges. We highlight two examples valid to an 
appropriate order in perturbation theory, uniform curvature gauge, which is as close to volume- 
preserving as possible, and a spatially-traceless uniform cold dark matter gauge which preserves the 
volume to linear order. We demonstrate the strong gauge- and frame-dependences in averaging. In 
traceless uniform CDM gauge the backreaction exhibits a strong ultra-violet divergence and can be 
tuned to an arbitrary magnitude with an appropriate choice of smoothing scale. In uniform curva- 
ture gauge we find that for a choice of Hubble rate locked to the spatial surface the backreaction 
vanishes identically, while for a Hubble rate defined from a fluid's expansion scalar the effective 
CNJ ■ energy density at the current epoch in an Einstein-de Sitter universe is O e ff ~ 5 x 10~ 4 , slightly 

bigger than but in broad agreement with previous results in conformal Newtonian gauge. 
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■ I. INTRODUCTION 



' q'| The gravitational field equations on cosmological scales are obtained by averaging the Einstein field equations of general 
q-i relativity. The effects of averaging (or backreaction) can have a significant dynamical effect on the evolution of the 
Universe and consequently on cosmological observations. The quantitative size of these effects, and their potential 
qualitative effect on cosmological observations, are currently of great interest (see, for example, [l|). However, in 
'— '■ order for the results of cosmological averaging to make any physical sense, it is absolutely crucial to have a rigorous 
| | (fully covariant) definition of the spacetime average of a tensor on a differential manifold. 

Cosmological perturbation theory provides a well-motivated paradigm in which to perform cosmological averaging. 
1 In perturbation theory we shall argue that a spacetime average is only well-defined when undertaken in so-called 
volume-preserving coordinates (VPC), or in the closely related comoving volume-preserving gauges (VPGs), which 
we define below. In fact, a VPG is ideally-suited to the study of spatial averaging in perturbation theory. 

We discuss two such gauges which are valid to an appropriate order in perturbation theory; the uniform curvature 
—i ■ gauge (which is uniquely well-defined in this respect and is as close to volume-preserving as possible) and the spatially- 
t-H \ traceless uniform cold dark matter gauge (which is adequate for calculations to second order in perturbation theory) . 
In fact, traceless gauges, including both of these, simplify the problem considerably since it reduces to an averaging 
of the product of linear perturbations and, at most, a second-order velocity. 
. An important measure of the effect of backreaction is the difference between the input Hubble rate and the averaged 
Hubble rate. We demonstrate that the cosmological backreaction from perturbations is strongly gauge-dependent, 
and clearly depends on the choice of frame for the averaged Hubble rate (e.g., the gravitational or the projected fluid 
frames). We shall show that in uniform curvature gauge the backreaction vanishes identically in the gravitational 
frame. We consequently argue that the definition of the Hubble rate should reference the fluid content of the universe 
and so it is necessary to use the projected fluid frame. We show that the backreaction in traceless uniform CDM gauge 
exhibits a strong dependence on choice of smoothing scale. In the projected fluid frame, the ultraviolet catastrophes 
in traceless uniform CDM gauge are exacerbated. 

We then turn to some quantitative numerical results. We find that in uniform curvature gauge the effective energy 
density of backreaction at the current epoch in an Einstcin-dc Sitter universe is of order 10~ 4 — 10~ 3 , and is slightly 
larger (by a factor of 2-5) than, but in broad agreement with, previous results obtained in conformal Newtonian gauge. 
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This slightly larger result perhaps suggests that backreaction within perurbation theory can be large enough to affect 
cosmological observations, but we believe the more important result is that this value is now on a significantly firmer 
basis than before. Of course, more significant quantitative results are possible in more realistic models in which the 
present-day universe is not well-described by perturbation theory. 

II. FORMALISM 

Three approaches to spatial averaging in perturbation theory and applicable in general coordinates have recently 
been presented • The relationships between these have been explored in 0, Q and are shown to tend towards 
agreement for extremely large domains. The explicitly gauge-invariant approach of [J, |5| has been further developed 
and applied to averages across the past light-cone in [JQ. However, our ultimate aim is to select a gauge in which 
the averaging procedure remains well-defined and so will focus on the formalisms of @> S @ ■ Likewise multiscale 
averaging models have been introduced (we see for instance [9l4l3|) but for simplicity we consider only single-scale 
models. 

A globally hyperbolic spacetime can be foliated with a family of 3-surfaces described by coordinates x^ = (t, x a ). 
The 4-vector normal to the surfaces is n M = (l/a)(l, — f3 l ) with normalisation n^n^ = —1 with the lapse a and the 
shift f3 l . The projection tensor onto the 3-surfaces, induced 3-metric and line element are then 

hfj,u = g^v + n^riy, hij = g^h^h^ ds 2 — (—a 2 + /3, ; /3 4 )eft 2 + 2(3idtdx l + h i} dx l dx° (1) 

where /3j = hij/3 J . The embedding of the hypersurfaces is described by the extrinsic curvature Kij = — (l/2)C n hij — 
— (1/2) (d/dt) — 2T)(i(3j). T>i is the covariant derivative on the 3-surface and brackets denote symmetrisation on the 
enclosed indices. 

Let I? be a finite domain lying on the inhomogeneous 3-surface, and let h be the determinant of the 3-metric. The 
volume of this domain (see, e.g. [1343) an d therefore the average over the volume of a scalar A can be defined by 

V v = [ Vhd 3 x, (A) = — I AVhd 3 x. (2) 

An average Hubble rate can be defined by 

3Hv - || - J (-2 1 '"' 1 -*) ^ x = (- aK + ( 3 ) 

with T> % = h^Vj. We dub this choice of Hubble rate the "gravitational frame" - it describes the Hubble rate defined 
with respect to observers comoving with the coordinate grid. 

It is possible instead to define an averaged Hubble rate from the expansion scalar of an observer, projected onto 
the spatial 3-surface, = h^ v u^ v . We term this the "projected fluid frame" and it was employed, for instance, in 
[3, @ . The Hubble rate is defined by 

m v = (a 2 e) . (4) 
Being defined from a physically-meaningful quantity this definition of the Hubble rate is perhaps to be preferred. 

III. AVERAGING IN AN ARBITRARY GAUGE 

The above equations are a rephrasing of the scalar 3+1 equations avera ged across a 3D domain and as such provide 
no new information. To make progress either some behaviour (as in [19(), or some underlying model (as in for 
instance @, [20l - [27j ) must be assumed. Cosmological perturbation theory provides a well- motivated toy case. This 
approach introduces a local scale factor a(t) and Hubble rate governing the unperturbed dynamics. One measure of 
the impact of "backreaction" is then to identify the backreaction as the difference between the input Hubble rate and 
the averaged Hubble rate. Of course, many other effects of backreaction will occur in cosmology In this paper, it 
is this Hubble- "backreaction" that we are interested in. 

The perturbed flat Robertson- Walker line element can be written 



ds 2 = a 2 (ji) (-(1 + 20)d? ? 2 + 2B i dr}dx i + (% + 2C ij ) dx l dx j ) 



(5) 
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where B. t = d L B - S, and C l3 = -ipSij + did 3 E + d^F^ + \h i:j with d l S l = d l Fj = 8%^ = S^h l3 = 0. The fluid 

velocity is Vi = diV + v\ . Perturbations are expanded to second order with </> = <j)(\\ + (1/2)0(2) and similar. We 
neglect everywhere products which are cubic or higher. 1 Indices are raised and lowered with the Kronecker delta. 
Derivatives with respect to conformal time will be denoted with an overdot and Ji = a/a is the conformal Hubble 
rate of the underlying model which obeys the Friedmann equations, 

H 2 = 8^ a 2£- l fl H = -^^(p !i) +3p (/) )+ia 2 A, (6) 
/ / 

Linear tensor modes are extremely small and decay with the expansion of the universe while, unless supported by 
an active source, vector modes rapidly decay. To a first approximation, scalar modes therefore dominate the linear 
perturbations [281 ] . At second-order all three types should in principle be considered, but it turns out that in the 
situations we study, the leading order contribution from vectors and tensors is quadratic, implying that to second-order 
in perturbations they can be neglected. 

The metric determinant of the perturbed model in a general gauge is 

Vh = 1 + C + ic* 2 - (-'(■,, . (7) 
If W(x) is a window function defining the domain, and 

Vjr = J T^(x)v//i (x)d 3 x = a 3 J W^(x)d 3 x (8) 
is the domain volume projected onto the FLRW background, the inhomogeneous domain volume is therefore 
V v = V T + a 3 J M/(x)C(x)d 3 x + a 3 J W(x) Qc 2 ( x ) - C iJ ( x )^j ( x )) d 3x 

with the integration taken across all x. When required to specify the window function we choose W(y) = exp(— y 2 ). 
Likewise, the average of a quantity A{x) — Ai(x) + A2(x)/2 is 

(AW) = ^J W(x)A(x) (1 + C(x)) d 3 x = ^ J W(x) ^x(x) + \a 2 {x) + ^(x)d(x)) d 3 x (9) 

while that of a product of linear perturbations Ai(x)Bi(x) is 

(A 1 (x)S 1 (x)> =^-J VF(x)A 1 (x)B 1 (x)d 3 x. (10) 

It is usual to then Taylor-expand the domain volume 

V v ' = V F l (l - ^ / ^(x)C(x)d 3 x + . . .) (11) 

where we truncate the expansion at linear order since corrections to averages of the type (A(x)) would enter at cubic 
order or above. 

Cosmological perturbations can be represented in Fourier space. Denoting an ensemble average with an overbar, 
the power spectrum of two linear perturbations is 

A^B*^) = ^(fc)A 1 (fc) J B*(fc)(2^) 3 <5(k - k') (12) 



with primordial power spectrum 



k sn.-l 



V(k)=A*[-) . (13) 



1 Note that this implies that if, for instance, the linear and second-order integrands are highly oscillatory, a dominant contribution could 
naturally arise at higher orders in perturbation theory. However, the study of such systems seems better suited to fully nonlinear 
rclativistic models. 
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The ensemble average of the domain volume is then 



Vv = Vr + a 3 J W(x)C(x)d3x+ V> J V{k) (\C{k)\ 2 - ^(^(fc)) y. (14) 

Likewise the ensemble average of a spatially-averaged product of linear perturbations becomes 

(^(x^x)) = - J V {k){A 1 (k)Bl(k)+c.c.)-. (15) 

The ensemble average of a linear perturbation is vanishing by definition. 
In the gravitational frame, the averaged Hubble rate ([3]) then becomes 

Hv = H + i ((7 - 2C* 'Cij) . (16) 

This corresponds with the large-scale limit of the the expansion scalar of the coordinate grid defined with respect to 
conformal time [29| , Q CO n( ~ 3 {% — ip — 2i/)ip J . The perturbed Hubble rate in the backreaction is a simple average of 
this quantity (denoted £ in Q). 



Squaring this Hubble rate gives 



and so the effective energy density 



Ul=U 2 + ^H(c-2C i W ij ) + ^(6f (17) 



H 2 n cS = ^H(C- 2C»C i0 ) + i (c) 2 . (18) 

It is important to note that in the gravitational frame the entire backreaction then depends solely on the choice of 
threading of the 3-surface, defined by the choice of Cij. The choice of slicing (i.e. of <f> and Bi) and the behaviour of 
any fluid content influences Ht> only indirectly through dynamics. 
In the projected fluid frame the averaged Hubble rate is instead 



H v = U + i (d i d l v + C + 2UC^ + -(j&diV + dtv&C - d t Bd l <p + ^Hd a Vd a V + 2C l3 d t d v - 2C l3 C tJ ) (19) 



and so depends both on the coordinates on the 3-surface, on the choice of slicing, and on the velocity of the fluid with 
respect to the background, v. The covariant velocity is V = v + B. 

Both these forms of the Hubble rate are trivially gauge-dependent. This is not surprising. The choice of gauge has 
governed our choice of 3-surface upon which to average and we should not expect gauge-invariance. Rather we should 
attempt to find in which gauges an average can be properly defined. 



IV. GAUGE CHOICES 



Coordinate freedoms allow us to eliminate two scalar degrees of freedom. Studies of backreaction have typically 
been in synchronous gauge with <f> = B = (e.g. [HI, [2 H3|) or conformal Newtonian gauge with B = E = 
(e.g-H, l24| - [27l [3TI [ipO . Uniform curvature gauge was proposed in [I?) but not examined in detail, although in a 
different formalism 33] a volume-preserving system with similarities to uniform curvature gauge was considered. 

In order to carry out any cosmological averaging, it is absolutely crucial to have a rigorous (fully covariant) definition 
of the spacetime average of a tensor on a differential manifold in order for the results to make any physical sense. 
There have been a number of recent approaches to this, including the exact macroscopic gravity (MG) approach which 
gives a prescription for the correlation functions that emerge in an averaging of the Einstein field equations [3^ - l4l1 |) 
(see also [l|). 

This is an absolute necessity for any results to be interpreted physically. In the context of this paper - averaging 
in perturbation theory - in general a spacetime average is only meaningful when undertaken in a volume-preserving 
coordinate system (VPC). 

In a VPC the volume of a domain is preserved as the system evolves in time. Closely related are what we term 
comoving volume-preserving gauges (VPGs) in which the volume of a 3-domain on an inhomogeneous surface evolves 
purely as a 3 (r/). In a VPC the time dependences cancel out when one takes an average - for the purposes of spatial 
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averaging, a VPG is then effectively a VPC. VPCs are employed in, for instance, Macroscopic Gravity [3J-|40( and a 
VPC at linear order was applied to cosmological averaging in 133] ■ VPCs are also explicitly utilized in an approach 
to averaging within unimodular gravity [42[. It is a central aim of this paper to study the significance of the VPGs, 
or whether it is adequate to average in another, practically convenient gauge. 

Written explicitly, the spatial average of a perturbation A(x) in arbitrary coordinates is given by 

Mfx s, = J ^(x)A(x)d 3 x + J VT(x)A(x)C(x)d 3 x 

[ [ '' /ly(x)d 3 x + /^y(x)C(x)d 3 x + /T^(x)(iC 2 (x)-Cy(x)C ^J (x))d 3 x■ 1 ' 

This takes on the simplest form when Cjj = 0. Tensor modes at linear-order are gauge- invariant and cannot be 
removed, but the gauge that imposes Cjj = (l/2)hij at an arbitrary order in perturbation theory is uniform curvature 
gauge. Neglecting tensor perturbations, 

W/ W( X )d 3 x , ( , 4(x)) = »^. 

This is then a VPG; the spatial surfaces align with the FLRW background, and the 3-volume expands only with the 
background and the only time-dependence in the average is that of the perturbation itself. This gauge is ideally-suited 
to the study of spatial averaging in perturbation theory. The tensor modes can be included. To second order, tensor 
and scalar contributions do not couple together, and the results of [13] can be directly employed: tensor modes from 
inflation will produce a baseline backreaction of order fi e ff « 10~ 14 . 

Uniform curvature gauge is the unique choice for a comoving VPG, but a less stringent alternative valid to second- 
order when averaging perturbations can also be found. If we assume that we can Taylor-expand Vj^ , we can reduce 
the spatial average to 

(a 3 / IU(x)A(x)d 3 x + a 3 / IU(x)A(x)C(x)d 3 x) U~^rJ C(x)d 3 x) 
(A(x)) — ^ ^ L (22) 

where we have truncated to second-order in perturbations. From this it is clear that if we choose a (spatially) traceless 
gauge, with C(x) = 0, the corrections to the volume are pushed to higher-orders in perturbation theory. Uniform 
curvature gauge is a trivial example of a traceless gauge, but a convenient, less stringent gauge can be found which is 
sufficiently close to a VPG and may provide a simpler basis for calculations than uniform curvature gauge. 
A traceless gauge is defined by 

C T = =J> 3ip T = d a d a E T . (23) 

The scalar transformation from an arbitrary gauge (written with a tilde) into this gauge is given by a generating 
vector £ M = (a, d l f3) with 

<t> T = 4> + a + Ha, C T = C + 3Ha + d a d a p = . . 

B T = B-a + P,5 T = 6-3H(l + w)a, v T =v-0. [ ' 

A choice containing no arbitrary constant is a uniform density gauge, where the spatial sections follow contours of 
constant density of some fluid. Aligning to CDM gives the generating vector 

a - , d a d a (3 = 3^ - ~S C + d a d a E. (25) 
Art 

In terms of Newtonian gauge quantities the traceless gauge perturbations are therefore 



5 c n 3 1/; 1,9 \H „ 

IPT = 4>N — ; E T = —j^ipr, 4>T = 4>N + ^7 I W + ~K z V cN - g^^cAT + cN 

Bt = v cN ~—6 cN , 5 a T = S aN - (1 + wa)8 c n, v aT = v aN - v cN 
on. 



(26) 



where we have used that in Newtonian gauge S c n = 3ipN + k 2 v C N. Note importantly that this gauge is therefore a 
uniform CDM density gauge comoving with CDM. 2 



2 It can be motivated as such; setting v c t = S C T = one recovers the above gauge with an arbitrary function of x appearing in Ex- Since 
this constant only influences the coordinates on the 3-surface and not the choice of slicing it can be fixed to ensure C(x) = 0. 
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It is interesting to note that the use of these gauges simplifies the evaluation of backreaction in an another important, 
practical manner: it removes inconvenient terms. Consider the averaged Hubble rate in the gravitational frame. This 
then contains the average of C. In an arbitrary gauge this contains spatial averages of both linear perturbations and 
second-order perturbations. While the ensemble averages of linear perturbations vanish, this still leaves an average 

across a second-order perturbation and products of spatial averages (Ci)- However, we only have a hrm method 

for evaluating the averages of products of linear perturbations. The second-order term is particularly problematic. 
Approximations for the nonlinear Bardeen potential in ACDM exist [f| [26|, |43[ and so, in principle, this term can 
be calculated in all gauges. These solutions are, however, valid only in matter-dominated universes (with or without 
a cosmological constant), and it would obviously be convenient to have a formalism that can be readily applied to 
universes with arbitrary fluid content. Choosing a traceless gauge immediately removes this term identically, and the 
form of the Hubble rate becomes significantly simpler. 

Similarly, in the projected fluid frame we must average across C(x) and d a d a v, and a traceless comoving gauge 
would appear to render the problem sufficiently straightforward. In uniform curvature gauge we still require knowledge 
of V2F- The desire, therefore, is that the almost-volume preserving traceless uniform CDM gauge can provide an 
"accurate-enough" approximation to the result in the genuinely volume-preserving uniform curvature gauge. 

In this manner one can motivate the choices of gauge we wish to study. Averaging in uniform curvature gauge 
is uniquely well-defined, and we consider it in both the gravitational and the projected fluid frames. However, the 
traceless uniform CDM gauge is also well-motivated and preserves the comoving volume up to linear order, and 
averages up to second order. In the gravitational frame, any traceless gauge - including both uniform curvature 
and traceless uniform CDM - reduces the problem to averaging the product of linear perturbations. In contrast, 
in the projected fluid frame a traceless gauge comoving with the fluid provides the most convenient choice. We 
will also consider Newtonian gauge in the gravitational frame to provide a direct comparison with previous results 

[I, US S0. 

As a final comment, we have been assuming that the corrections to Vx> are small, and that we can take a Taylor 
expansion of the inverse volume. Physically one would not expect an integral across a gravitational potential to give 
a large result, but the correction integrals in the volume are frequently divergent in the infra-red, and potentially also 
in the ultra-violet (see for instance the power spectra in [45l - l47l | ; calculations of the volume reduce to integrals across 
such spectra). These divergences contribute to the spatial averages, and could render the Taylor series physically 
valid, but mathematically suspicious. Since it is physically implausible that metric perturbations cause the volume of 
a domain to become arbitrarily large we follow standard procedure and neglect this issue. Nevertheless, it is worth 
noting that the only gauge in which this does not arise is the uniform curvature gauge in which Vp — Vjr 1 neglecting 
tensor modes. Preserving the consistency of the formalism may then demand that we work in this gauge. 



V. DYNAMICS 

Up to this point our treatment has been applicable to any perturbed flat FLRW universe. For simplicity 3 we choose 
henceforward to work in a pure dust Einstein-de Sitter universe with £l m = 1, fl\ — and % — 2/77, unless stated 
otherwise. When necessary we choose fib — 0.05 and h = 0.704. The amplitude of primordial perturbations will be 
A+ = 2.42 x 10~ 9 , and we employ a Harrison-Zel'dovich spectrum with n s = 1 for simplicity. A small spectral tilt 
will not significantly change our answers. 

This model serves as a reasonable approximation up until recent redshifts. Analytic solutions for the Newtonian 
potentials in a dust-dominated universe (with or without a cosmological constant) at both linear and second-order 
were derived in [43| and applied to the averaging problem in 0, [26jj, and we quote the results directly here. 

In Newtonian gauge the linear potentials, fluid density contrast and velocity are given by 

4>in = 4>in = .9(»?)0o(x), S 1N = -2(f> 1N + ^r] 2 d a d a <f>iN -wPin, v 1N = -\r)(j) 1N - \if<j)i N (27) 

where g(rj) is the growth function and 0o the value of the Newtonian potential at the present epoch. At second-order, 
the potentials are given by 

i>2N = A 1 (r))<f>l + A 2 {r } )x2((l>o) + A 3 (ri)x3{<t>o) + Mv)d i <f>odi(t>o, (28) 

(f>2N = Mv^l + Mv^iM + MvhsiM + MvWfodifo (29) 



And to ease direct comparison with previous work. 
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with (j>2N 7^ 1P2N due to an effective anisotropic stress arising from products of linear perturbations. Here A n (rj) and 
A n (Tl) are functions of time related to the expansion of the background which can be found in [43j . while Xniv) arc 
second-order products of gradients and inverse Laplacians of <po, 

5-2 /oii Q J. \ oQ-4a flj /au a i. \ /j, \ _ a— 2/ 



X2(0o) = (SVo^o) - 3d" 4 ^ {d*<t> dj<t> ) , xsW-o) = cT 2 ^ (Pfadjfo) 



(30) 



To find the velocity in uniform curvature gauge we also require the Laplacian of the scalar velocity, which can be 
found from the momentum constraint [48| , 

^Ga 2 {p + p)d l d t v 2N = -d l d^2N ~ Ud l d^ 2N + 8irGa 2 (p + p)d l {(pi N diV 1N ) + 16nGa 2 d l {^ 1N d l v 1N ) (31) 

-8nGa 2 d l [(Spur + 5 PlN )d lVlN ] - AO 1 {^in^in)' + 2d 

In EdS g(rj) — 1 and the solutions reduce to 

4>1N 



di4>\N 



IN = (j) , S 1N = -2(f>l + ^rj 2 d a d a <f)o, v 1N = ~T]<j)o, 



ip 2 N = -205 - -X2(M + B 3 (rj) ( X3(M - j^d^odtfo 
"2\ '^"n + -V2("i.i + B-[ij) I x 3 (<M - ^fodifo 



1p2N = 4>: 



'2N 



with 



i (n. 



14 



V 



1 / Vm 
7 \ T, 



5 2 

2l"' 



(32) 
(33) 

(34) 



Here rj m is an early time deep in matter domination at which the Newtonian potentials are initialised and we focus 
on the regime 77 3> rj m . </>o can be readily recovered from a Boltzmann code. The divergence of the velocity potential 
is then 



-d a d a v 2N = ^V 2 d% (d^odacbo) - ^V 2 9% (difa&fa) - ^d^od^o + ^ 2 d l {d i( j> d a B M + 2 X3 (cf> )- 
The linear gauge transformation from Newtonian to uniform curvature gauge is generated by the 4- vector 



1plN_ 

H • 



2 #o,0 



With this transformation vector the linear uniform curvature quantities are readily found to be 



'IF = 2^0, 



"^#0, SlF 



1 



5 + ~k 2 r) 2 ) 0o, vif = -~#o, Vif = -~#o- 
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G 



The gauge transformation for the second-order velocity potential can be written [29j as 

d a d a v 2F = d a d a v 2N - d a d a p 2 + d k x v k 

with 

X\ = -2 ai di (v 1N + Hv 1N ) = # o 9i0o, h = --^d^x 1 ' + \d- 2 xt 

and the gauge function Xij is 

2 / ■ \ 2 1 

Xij = -jji>iN [tpm + 2UiPin) Sij - -^di^iNdjiiiN = -4i/; 1N Sij - -rfditpindjiiiN . 



H 



(35) 
(36) 

(37) 

(38) 
(39) 

(40) 



After some manipulation this gives the uniform curvature gauge velocity in a pure dust universe as 

d a 8 a v 2F = (d^odafo) - ^fd l dj (difo&fa) + \^d l (drf Q d a d a fo) 

5 17 

+ 4?7X3(0o) + Tl<j)od a d a (j)o - —rid 1 <M*<?V (41) 

Finally, we need the curvature perturbation, lapse and shift in traceless uniform CDM gauge. From equations (|26[) 
these are 



1>rr = fto - ^v 2 d a dM, <f> 1T = ^rfd a d a cf> Q , Bit = ~^V 3 d a d a <f> . 
3 18 36 36 



(42) 
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VI. THE GRAVITATIONAL FRAME 
A. Uniform Curvature Gauge 

A first step towards the use of uniform curvature gauge in this frame was presented in Q. However, in that study 
the authors only presented general forms and did not interpret the results. From equation (|16p . the averaged Hubble 
rate and effective energy density in uniform curvature gauge are then 

H v = H, H 2 n cS = 0. (43) 

The backreaction in uniform curvature gauge is identically zero, to an arbitrary order in perturbation theory! This 
result is contrary to the claim in Q that a gauge cannot be found that removes the backreaction. 4 

Presented in this manner, this result is trivial: the impact of scalar and vector perturbations on a Hubble rate in the 
gravitational frame vanishes identically, because we are working in a gauge with vanishing spatial scalar and vector 
perturbations. It is straightforward to interpret this result: the averaged Hubble rate in the gravitational frame is 
defined by the change in the volume of the domain. Since we are in a comoving volume-preserving gauge, the volume 
expands only with a 3 , and so the averaged Hubble rate is given purely by the input Hubble rate. 

This vanishing answer contradicts expectation - while certainly one might argue the backreaction (that is, the 
"Hubblc-backreaction" ) from perturbations should be small it cannot be expected to be identically zero at second- 
or higher-orders in perturbation theory. Given that uniform curvature gauge provides the best-motivated system in 
which to average, this suggests that the gravitational frame is ill-suited to studies of backreaction. Certainly it is 
difficult to connect the Hubble rate averaged in this frame with any physical quantity. 



B. Traceless Uniform CDM Gauge 

The above conclusions follow trivially from the definition of the Hubble rate in uniform curvature gauge, but as argued 
earlier we can relax the gauge constraint slightly and employ a traceless gauge. In the traceless uniform CDM gauge, 



4 In that study the authors did not address the effective energy density and instead demanded that each individual backreaction term in 
the "Buchert" approach vanish, and a gauge cannot be found in which that is the case. A gauge can, however, be found in which the 
combination of the backreaction terms vanishes, which we explicitly demonstrate in Appendix fX] 
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the averaged Hubble rate and effective energy density become 



H v =n- -(cVCy), H 2 n cS = -){((■■>(;,). (44) 

Since Cy- = —ipT&ij + didjEx and Sipr = d a d a Ex, the ensemble averages can be rewritten using equations (|15l) to 
become 

^ = n ~ i"/ k2v{k) l<Mfc)|2 f - ir? 3 / kiv{k) ^ o(fc)|2 T- (45) 



We have written the curvature perturbation in this gauge in terms of the Newtonian gauge quantities using the gauge 
transform l[4"2"]). 

We recover <fio from a modified version of the CMBFast code [49| itself based on Cosmics [5(|, but we can gain 
insight examining the zero baryon transfer function found in [5lj , 

L 731 

T = 2 , g = fc Mpc /i-'Gly/r, r = ^ = /», Lo=M2e + 1.89), C = 14.2 + — — — . (47) 

This provides a good approximation to the numeric 4>q; with fib — 0.05 the baryon oscillations and small-scale 
damping from the baryons is relatively minor. The analytic form shows that there is an ultra-violet divergence 
in the term proportional to k 4 ; on large scales the integrand scales as ~ (\nk) 2 /k which produces a logarithmic 
divergence ~ (ln/c) 3 /3. We control this as in (26|, smoothing the gravitational potential in real space with the 
window function W(x/Rs)- The smoothing scale Rs is arbitrary. A well-motivated choice would be the Silk scale, 
Rs = (fcsak) 1 ~ 6Mpc where the numerical value is for an Einstein-de Sitter universe. 

Figure [T] shows the integrands, generated by combinations of the Newtonian gauge gravitational potential, which 
dominates on large scales, and the Newtonian gauge density contrast, which dominates on smaller scales. It is 
clear that the integral does not contain an infra-red divergence, and the integrand with and without the small-scale 
smoothing is plotted, taking Rg = 6Mpc. 

We solve the integrals numerically. For R$ = i?suk we have 

Uv ~ U = -0.409, Q eS = -0.818 (48) 
H 

(Due to the additional small-scale damping in the baryon case the results found from the Eisenstein and Hu zero- 
baryon transfer function (|47|) are slightly larger: (Hx> — H)/H = —0.511, fi e ff = —1-02.) The unphysically large size 
of these results strongly suggests that in this gauge the perturbations must be controlled on a larger scale. 

The left panel of Figure [5] shows the fractional shift (Hx> — T-Cj/H as a function of smoothing scale Rs for the test 
Einstein-de Sitter universe. The dependence is extremely strong. For Rs < 4Mpc the fractional change to the Hubble 
rate is significantly larger than unity. The magnitude of the correction decays monotonically as Rs — > oo; since the 
entire integral is smoothed there is no asymptote. It is important to note that while there is a dependence on the 
smoothing scale Rs, there is no dependence on the averaging scale 1Zt>- 

In the right panel of Figure [2] we plot instead O e ff = {R-t> — 'H 2 )/'H 2 . Qualitatively the behaviour with Rs is the 
same as that for the Hubble rate. Both (Hv — R)/T-L and f2 e ff are negative in this gauge and frame. To obtain 
n cB < 10 -2 would require R s > 24Mpc. 

The arbitrariness of the results is extremely unsatisfying and stems from the strong ultra-violet divergences. The 
automatic conclusion is that contrary to expectation this gauge in the gravitational frame is not well-suited to 
calculations of backreaction. Moreover, since O e g = in the uniform curvature gauge, traceless uniform CDM gauge 
certainly can not be used to approximate the more well-defined choice of gauge. However, it is interesting to note 
that were we to take an infra-red divergence that appears in the domain volume seriously, and smooth the ultra-violet 
divergence in the integrals above, then the effective energy density is driven to zero. 
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101) 

Rs (Mpc) 



FIG. 2: Modification to Hubble rate (left) and effective energy density (right) at the current epoch as a function of smoothing 
scale Rs- The domain radius is Rd = r/o/2. 



C. Conformal Newtonian Gauge 



The final gauge we consider in the gravitational frame is conformal Newtonian gauge. With C'ij = —ipN^ij, the 
average Hubble rate and effective energy density are 



1 



1 • 



Hv = H + -( IplN + ~1p2N - 6tplNlplN ) = -(i>2N), 



1/7" 



n 2 n 



cir 



(49) 
(50) 



where we have used that ip\N — to simplify the forms considerably. Using the analytic expression for ip2N m equation 
reduces this to 



^n) = ^B 3 (V) J k 2 V(k) \Mk)\ 2 y « ^/ *? V W l^°( fc )| 2 T 
where we have used B3 (10/21)^. We therefore have the average Hubble rate and effective energy density 

n^ = n + ^ v f k 2 r(k) \Mk)\ 2 y , n 2 n cS = A f k 2 v(k) \Mk)\ 2 y . 



(51) 



(52) 



These integrands are plotted in Figure [T] and contain neither infra-red nor ultra-violet divergences. Since the term 
of the form (•) (•) is vanishing in matter domination there is no dependence on the averaging domain; in a ACDM 
universe, however, this term would be non- vanishing and there would be a weak dependence on 1Zt>- 5 

It is also clear that in matter domination the effective energy density of the backreaction becomes a constant, in 
qualitative agreement with 0, [25|, [26[ . 



5 Strictly speaking this contradicts a statement in Q that the backreaction in gravitational frame does not depend on the choice of 
averaging domain. However, even in a ACDM universe ip <C if) and the contribution will be entirely negligible, and for all practical 
purposes the conclusions of that paper are unchanged. 
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From the form of the integrands the modifications to the Hubble rate in Newtonian gauge will be both significantly 
smaller than, and of an opposite sign to, those in uniform traceless CDM gauge. For the test EdS cosmology, we find 



Hv-H 



3.44 x 1CT 7 , fieff w 6.87 x 1CT 7 . (53) 



H 

Using the zero-baryon transfer functions yields the slightly larger (Hv — H)/'H = 4.4 x 10~ 7 and f2 c g «9x 10~ 7 . 

VII. THE PROJECTED FLUID FRAME 

Averaging in the projected fluid frame tangles together fluid and metric quantities and it is not possible to significantly 

simplify expressions employing the dynamical constraints. In particular, the term Ht>,i — (d a d a v + C + 21-LC^j is 

present at both linear and second-order in perturbations and does not readily simplify. For instance, at linear order 
the Laplacian of the velocity can be replaced with a combination of the lapse, shift and curvature using the momentum 
constraint ip + H(f> = — 47rGa 2 (u + B); removing the velocity then results in the average of a combination of <j>, B, ip 
and E, which will not cancel. The situation at second-order is significantly more complicated. 

Uniform curvature gauge remains the uniquely well-defined gauge in which to perform a spatial average. It can 
then be argued that a calculation of the backreaction in this frame and gauge is the best we can hope for within the 
confines of cosmological perturbation theory. 

In the gravitational frame, we found that the traceless uniform CDM gauge is not a good approximator for the 
uniform curvature gauge, despite being almost volume-preserving. In the fluid frame there is an extra motivation 
for studying this gauge: we can calculate the backreaction true to second order using only linear perturbations since 
vt = Ct = 0. If the uniform traceless CDM gauge is a reasonable approximation to the uniform curvature gauge, it 
is then significantly more straightforward to find results in this gauge. 

A. Uniform Curvature Gauge 

In uniform curvature gauge, the averaged Hubble rate (|19l) simplifies slightly to become 

Hv = H + i (d a d a v 1F ) + i (d a d a v 1F ) + i (to d a d a v 1F - d a tod a B 1F + ^Hd a V 1F d a V 1F ^ . (54) 
The ensemble average of the Hubble rate and the effective energy density are therefore 



Hv = H + Ud a d a v 1F + ^d a d a v 2F \ + Utod a d a v 1F - d a tod a B 1F + ^Hd a V 1F d a V 1F ). (",) 



n 2 n cif = ^n/d a d a v 1F + ^d a d a v 2F 



+-H(tod»d a v 1F - d a tod a B 1F + ~Hd a V XF d a V 1F j + -(d a d a v 1F y. (56) 



We consider these averages term by term. Since the ensemble average of a linear perturbation vanishes, {d a d a v\ F ) = 0. 
From equation (|4"Tj) we can also see that 



( ^vxsiM - ??0o9 a a a 0o - jr; '/<>'•■ ' ) (r>si 



The first of these terms is 



/io ,„.„ , A 10 
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Transferring this to Fourier space reveals that this term vanishes on ensemble averaging. The second term is 



1 „„.„ ,„ , „ , A 1 



which is readily seen to also vanish on ensemble averaging. The third term, 



1 „ , „ „ , A 1 



gf? 3 ^ (O^od-daM ) = -^(d^od-da^o + d^d-dafo) (61) 



also vanishes. In 2_6j it is shown that 



and since 



we can see that 



(Xs(0o)) = gWoft^o) (62) 



dk 

(d»<hd a <M = -(fod a d a cbo) = I k 2 T(k) \Mk)\ 2 y (63) 



{d a d a v 2F ) = 0. (64) 
The ensemble averaged Hubble rate then reduces to 



25 f,J fl mu „M2d* 



Hv=H + - V k 2 V(k) \Mk)\ y (65) 



In [2g| it was additionally shown that 



dfc 



{d a d a Af = I fc 4 |^L(fc)| 2 W 2 {kR v )^- (66) 

where Rx> is a length scale characterising the averaging domain. Using this we can see that the effective energy density 
in uniform curvature gauge and the projected fluid frame is 

u 2 n cS = y J k 2 v(k) \Mk)\ 2 Y + h v " J k4v{k) l0o(fc)|2 w2 ( kR v)jr- (67) 

The only scale dependence in the solution enters in this final term - for a large enough volume, the integral is driven 
to zero and the effective energy density is governed by the first term. 

The right panel of Figure [T] shows the integrands of H-v and of f2 e g. We can see that the impacts on the Hubble 
rate and its square are both positive, are cleanly under control, and will be smaller than those in traceless uniform 
CDM gauge. Setting the domain scale to the Hubble scale gives 



U-d-U 
H 



1.8 x 1CT 4 , Q eS = 3.61 x 1CT 4 . (68) 



Using the zero-baryon transfer function yields the slightly larger (H-p — H)/H — 2.3 x 10~ 4 , il c g = 4.6 x 1CP 4 . This 
result is directly comparable with that in Q which performed the equivalent calculation in conformal Newtonian 
gauge, with the result ft e g » 4 x 10~ 4 , (H-p — T~L)/T~L ~ 2 x 10 -4 for an Einstein-de Sitter cosmology. It is also 
approximately in line with earlier calculations such as those in [24], [H, [52[ • 

Figure [3] shows the dependence of the effective energy density on the averaging domain scale. At Rp = 16Mpc the 
effective energy density is O e gf = 8x 1(U 3 , decaying to fi e ff = 4.4 x 10 -4 at Rp = 64Mpc. In this gauge it is then 
possible to identify a loose "homogeneity scale" at Rp, ~ 150 — 250Mpc, above which the backreaction becomes scale- 
independent. This agrees well with the calculation in [(| [2(|, which identified a similar scale in conform Newtonian 
gauge. 
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FIG. 3: Effective energy density of the backreaction in the projected fluid frame and uniform curvature gauge, as a function of 
domain radius Ri> from accurate numerical calculations (black) and the zero-baryon approximation (blue). 



B. Traceless Uniform CDM Gauge 

The traceless uniform CDM gauge comoves with the CDM, and is volume-preserving to first-order in perturbations. 
Selecting this gauge, the averaged Hubble rate is 

Uv = H + ~ (hid a Bd a B - d a cf>d a B^ - | (d ab C ab ) (69) 

and we can therefore calculate the result up to second-order in perturbation theory employing only linear perturbations. 
The second term is identical to that in the gravitational frame, and so 

Hv=n + n Giav . + i (^W a Bd a B - d a cl>d a Bj . (70) 

Using the perturbations given in (|4"2j) we find that 



-Ud-Bd a B - d a (f)d a B J = ^ J k b \M —. (71) 

This contains a severe ultra-violet divergence which we control with a smoothing scale Rs- The severity of the 
divergence suggests we must smooth on much larger scales to control otherwise divergent results. 
The averaged Hubble rate and effective energy density are then 

H^-n - -^if J ' k e V(k)\cb Q \ 2 W 2 (kR s )^ -^Grav., (72) 

n 2 n cB = [ k 6 v{k)\<f> \ 2 w 2 (kRs)^--H 2 n eH , Glav .. (73) 

The integrands are plotted in Figure [T] On super horizon scales these agree with the traceless uniform CDM results in 
the gravitational frame, but on subhorizon scales the ultra-violet divergence is extremely notable. Even smoothed at 
Rs = 6Mpc, it is clear that the results in the projected fluid frame will be orders of magnitude larger than those in 
the gravitational frame. There is also a sign change, suggesting that the sign of the backreaction, at least, will be in 
agreement with the more controlled calculations in uniform curvature and conformal Newtonian gauges. Evaluating 
the present-day backreaction at Rs = 6Mpc, 

nv ~ n = 3.9 x 10 4 , n cS = 7.82 x 10 4 ! (74) 
rL 
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The solutions as a function of Rs, for Rx> — T)q/2, are presented in Figure[2] As might be expected, for a sufficiently 
large Rs - a smoothing scale approaching that of the Hubble scale itself, Rs > 2000Mpc - the results in the two 
frames coincide with one-another. On smaller scales the effective energy density is indeed positive, but we require 
a smoothing scale Rs > 60Mpc if we want il c s < 1 - to ensure fi e ff ^ 0.01 we need Rs > 130Mpc! The extreme 
divergence, and the smoothing scales required to control it to recover meaningful results, suggest that the 3-surface 
and its perturbations in this gauge are badly-suited to the study of backreaction. It is certainly possible to argue that 
a comoving uniform density gauge is not well-adapted to the recent universe; the recent universe contains both large 
velocities and high densities, implying that the hypersurface and its embedding grow increasingly contorted. In any 
event it is certain that we cannot use the traceless uniform CDM gauge as a convenient substitute for the uniform 
curvature gauge. While in an EdS universe a choice of Rs ~ 243Mpc recovers equivalent results, we can not expect 
the same to hold true in more realistic universes - the smoothing scale required for consistency will not remain the 
same. 



VIII. DISCUSSION 

In this paper we have shown that the cosmological backreaction is both highly gauge-dependent and highly frame- 
dependent. We have advocated the use of a (comoving) volume-preserving coordinate system as the system in which a 
spacetime average is well-defined, and argued that given the smallness of the tensor perturbations uniform curvature 
gauge provides a surface to average across which preserves the comoving volume, and in which the time-dependences 
in the averages cancel. 

We then motivated an alternative choice of a comoving VPG, a gauge chosen to comove with surfaces of uniform 
cold dark matter density, with the spatial coordinates chosen such that the trace of the spatial metric vanishes. While 
this gauge does not preserve the comoving volume at second order in perturbation theory it does preserve it at linear 
order, which is adequate for calculations to second order in perturbation theory. A convenient feature of this gauge 
is that one can solve the backreaction to second order employing only linear perturbation theory. 6 

We compare both of these gauges against the conformal Newtonian gauge, employed in |jj [H, [2(| HH, [32J. In 
the gravitational frame, the backreaction in uniform curvature gauge vanishes identically. Since it is in the uniform 
curvature gauge that the backreaction should be defined, this suggests that the definition of the Hubble rate from 
the expansion of a 3-volume is too restrictive; with no reference to the fluid content of the universe it is also hard to 
recover meaning from the results. Further, while the backreaction in traceless uniform CDM gauge exhibits a strong 
dependence on choice of smoothing scale Rs , results in the gravitational frame do not exhibit any dependence on the 
choice of averaging scale Rx> which, as pointed out in [6|, [26[, is rather unnatural. The effective energy densities in the 
traceless uniform CDM and conformal Newtonian gauges possess different signs, and to ensure they are of equivalent 
size we must use a smoothing of the order of hundreds of megaparsecs. It is clear that it is not possible to use the 
traceless uniform CDM gauge as a simple alternative to uniform curvature. 

In the projected fluid frame, the ultraviolet catastrophes in traceless uniform CDM gauge are exacerbated. To 
ensure rt c g < 1, as seems reasonable, we must smooth perturbations on scales smaller than Rs = 64Mpc! In contrast, 
in uniform curvature gauge we do not need to smooth perturbations on small scales, and we find O e ff ~ 4 x 1CT 4 . This 
is in line with the results of [f| [26| in conformal Newtonian gauge, and also agrees with previous order-of-magnitude 
estimates such as in [25| . The use of uniform curvature gauge does, however, require knowledge of fluid velocities at 
second-order in perturbation theory, which is in general non-trivial. 

The effective energy density in the projected fluid frame the backreaction also exhibits the expected dependence on 
the averaging scale Rt>, with the impact decaying and asymptoting to a constant as Rx> — > Hq 1 ■ This is in agreement 
with the behaviour noted in [f| [26[ . In these papers it was also stated that the two frames should agree when the 
averaging scale is on the order of the Hubble scale. We confirm this for traceless uniform CDM gauge, but only if the 
ultra-violet smoothing scale is itself approaching the order of the Hubble scale. 

In uniform curvature gauge, however, we find that the effective energy density tends towards a a constant, finite 
value in the projected fluid frame as the averaging radius grows to infinity, which contrasts with the identically 
vanishing result in the gravitational frame. In the gauge in which the averaging is properly defined, the results in the 
two gauges will never coincide no matter how large the averaging domain. This forces us to choose a frame in which 
to work; since it is defined from physical quantities that have meaning for an observer we advocate the use of the 
projected fluid frame. 



Note, however, that this relies on it being possible to fix the gauge at second order in the same manner as we have at first order. While 
this seems likely, it has not been demonstrated. 
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We have presented a calculation of the cosmological backreaction in pure matter universes in the uniform curvature 
gauge, in which averaging is well-defined. The effective energy density of backreaction in this gauge agrees well 
with previous calculations in conformal Newtonian gauge. An alternative, which is well-defined to second-order in 
perturbations, does not provide consistent results. This gauge also ensures that corrections to the 3-volume remain 
formally small and that Taylor expansions remain valid, which is not the case in either of the alternative gauges. While 
the backreaction remains of order 10~ 4 — 10~ 3 , as in previous calculations, this value is now on a significantly firmer 
basis than before. The result is also slightly larger than the estimates in, for instance, 0, [25[ and the calculations 
in [(| [2^, While the direct impact from perturbation theory is still relatively minor, it is not so clear that it 
can simply be neglected. The present day universe is not well described by second-order perturbation theory. If 
perturbations induce backreactions at the order of 10~ 4 — 1CP 3 , larger inhomogeneous structures could be expected to 
have a larger impact - conceivabl y o f the order of > 10 -2 and equivalent to the energy density of baryons themselves 
(see for instance 0, [H,[5j] and [22j]). 

The result in a A or 0CDM universe will be somewhat less due to the washing out of structure from dark energy. 
Our result is then an upper limit on the present-day impact of second-order perturbations on the background. A 
more comprehensive study would require examination of the deceleration parameter or on other measures (such as 
the variance of the Hubble rate, considered for instance in [6, 55, 56]). Further progress, valid in the present universe, 
will then likely require the study of fully non-linear solutions to GR. 
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Appendix A: The Backreaction in the Buchert Approach 

In this appendix we prove in the context of the standard "Buchert" approach that the gravitational frame back- 
reaction vanishes identically in uniform curvature gauge. It is standard to connect the averaged Hubble rate to the 
fluid content of the universe. Applying the averaging procedure to the Hamiltonian constraint and evolution of the 
extrinsic curvature produces Friedmann-like (so-called "Buchert") equations in the domain 

Hi = ^ > Un, lS) ) - - > A - - | R P Q' p - (, > X< ' 1 . (All 



3 

a-D 3 



£ (a 2 ( PU) + 3^)) + | (o?) A+Uvl + Qi-Sj:^). (A2) 

/ V / 



Qp, V^p , TZ-p and Tx> are, respectively, the kinematic backreaction, dynamic backreaction, averaged curvature and 
fluid tilt, which corrects the fluid quantities between the surface orthogonal to n M and the rest-frame of a fluid with 
4- velocity u^y and are given in [3]. is the rest-frame density of a fluid and its rest-frame pressure. The 
effective energy density is then 



H 2 n cS = - | *> Wp' Rp -Q'p] ■ l A3) 

Using Bi = diB, u5> = 9 l i>( Q ), the correction terms become 

n v = 0, Q v = -AU (d%B) + ((d l diB)(d J djB) - {d%B){didi B)) , 
T U) = ^a 2 p (/) (S U) + 20) + ^a 2 p (/) (2# (/) + d^B + (1 + w^V^dW^) (A4) 

C v = \a 2 k (cf>) + ra 2 A (d l Bd t B) . 
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Expanding the perturbations into first- and second-order components, the effective energy density (|A3I) becomes 
^a 2 Pcfr = ^ P(f) + 2 ^(D + \ 6 im + ^(2)) + \a*A ( 2 ^(D + 0(2)) + (2B (1) + S (a) ) ^ 

+ ( ^T^XX) ( 2 ^(i)<W) + ^(i)^(i) + (1 + «J (a) )ftV( , 1) a i y (a)1) ) + ia 2 A9 l i? (1) 9,i? (1) ^ 

- ((3*^(1)) (^d)) - (^B(i)))) • (A5) 

The second-order perturbed Hamiltonian constraint in uniform curvature gauge [48j is 

ZUd%B (2) +m 2 4> {2) +^Ga 2 Y J 7{j)hf,2) = 

/ 

+6H 2 (# 2 1} - tfB^diB^ + m^x^diB^. (A6) 

Employing this to eliminate d l diB(2) in the effective energy density, and using the Friedmann equation to absorb 
terms proportional to A leads ultimately to 

H 2 n cS = l^a 2 Y^-p {f) 5 {L1) +2H 2 (1) + ^Hd%B {l ^j (A7) 

The above form of the effective energy density contains only first-order perturbations, for which we have a complete 
and straightforward theory. In particular, we have the Hamiltonian and momentum constraints 

3H 2 <j> (1) + nd%B w = -47rGa 2 ^7%)%i), Ucf> = -AirGa 2 £(1 + w (f) )p U) V (M (A8) 

/ / 

which can be used to easily eliminate the fluid quantities in terms of metric quantities. Use of the Hamiltonian 
constraint removes the average of first-order perturbations that appears in T-L 2 fl c ff, which becomes 

n 2 n cS = ^ a 2 ^p (/) (20 (1) % 4) + (i + W(/) )5 l i? (1) a l % 1) ) 

+m 2 4> 2 {1) + ^n4> W d%B (l) + ^Hd^diB^J . (A9) 

The Hamiltonian and momentum constraints give 

20(i)^a 2 ^P (/ )%,i)+4H 2 2 1) + ^ (1) 9 l 9 J S (1) = 0, 
6 f 6 

0*B (1) ^a a V(l + w s )-p {J) d % V UA) + hid'B^dicP^ = 0. (A10) 
6 f 6 

A quick examination of the effective energy density of the backreaction quickly confirms that 

H 2 n c{i = 0. (All) 

We have verified, within the full Buchert approach, that the gravitational-frame backreaction up to second-order in 
perturbation theory in uniform curvature gauge vanishes identically! Note that in forming this conclusion we have not 
transferred the system into Fourier space nor have we taken an ensemble average - the conclusion follows inevitably, 
in real space, for scales on which second-order perturbation theory is valid, and for any admixture of fluids. 



